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The 0{N) invariant quartic anharmonic oscillator is shown to be exactly solvable if the interaction 
parameter satisfies special conditions. The problem is directly related to that of a quantum double 
I ^1, well anharmonic oscillator in an external field. A finite dimensional matrix equation for the problem 

f-^ is constructed explicitly, along with analytical expressions for some excited states in the system. The 

corresponding Niven equations for determining the polynomial solutions for the problem are given. 

a 

i—j , Keywords: Quartic anharmonic oscillator; quasi-exact solutions; Niven equation 

PACS numbers: 03.65.Ge, 02.30.Mv, 02.20.Qs, 02.60.Lj 

(N 
> 

Anharmonic oscillators continue to be of interest because of their enormous phenomenological significance. '^1 In 
(T^ ] particular, the quartic case has been studied extensively in different contexts for classical and quantum mechanical 
^—\ ' systems, especially in one- and two-dimensional spaces. '^"^l This is of particular interest as it can help shed light on 
S-J ' the structure of the corresponding t/i^ quantum field theory. Up until now, most of the results have been based on 
CsJ ' a strong coupling expansion, [*1 variationally modified perturbation methods, 1^1 or lattice methods, [^°1. In each case 
only numerical results are accessible. 

Analytical solutions of the anharmonic oscillator problem for polynomial potentials of degree 4fc + 2, where fc is a 
positive integer, and some model specific restrictions have been reported. For example, the sextic and dccatic cases 
d i have been studied rather extensively.'"'^^"^^' The k = 1 case was first studied in [11,12]. This class of anharmonic oscil- 
lators is conditionally solvable in a coordinate representation, in which finite polynomials multiplied by an exponential 
can be adopted as the ansatz for wavefunctions. The quasi-exactly solvability of these cases is explained in detail in 
[23]. However, the important quartic case, together with polynomials of degree 4fc for fc > 1 is absent from this list of 
quasi-exactly solvable cases on the full line. The quartic anharmonic oscillator in one dimension on the half line with 
potential V{x) = ax^ -\- fix"^ and x > was first studied by Flessas et al in [24, 25] in terms of an integral expres- 
sion. It was commonly believed that the lowest-degree one-dimensional quasi-exactly solvable polynomial potential 
is sextic. This belief was based on the assumption that the Hamiltonian must be Hermitian. However, it has been 
discovered recently that there are huge classes of non-Hermitian, parity and time reversal invariant (PT-symmetric) 
Hamiltonians whose spectra are real, discrete, and bounded below. Replacing hermiticity by the weaker condition of 
PT symmetry allows for new kinds of quasi-exactly solvable theories. '^^^^^1 In [26, 27], a new two-parameter family 
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of quasi-cxactly solvable quartic polynomial potentials with V{x) = —x + 2[ax^ + [a — 2b)x + 2i{ab — J)x, where 
a 7^ 0, 5 7^ 0, and J is an integer, is introduced. In [28], the real spectrum of the potential V{x) = x'^{ixY with v > 
was studied. In this paper, we show that exact solutions of the 0{N) invariant quartic anharnionic oscillator in the 
spherical coordinates can be exactly obtained when the interaction parameter satisfies some special conditions with 
an ansatz similar to that used for the scxtic casc.I^'^1 

In A^-dimensional space, the 0{N) invariant Hamiltonian for the quartic anharmonic oscillator with unit mass can 
be written as 

i 

where V{r) is a polynomial in r = (^^ xf) of degree 4. Because of the underlying 0{N) symmetry, the Schrodinger 
equation of (1) can be transformed into the following form in A^-dimensional spherical coordinates: 



-iA + y(r))*(x,) = S*(a;.) (2) 



with *(.Ti) = Ri{r)Hi{n) and 



, , , / (f N-ld l(l + N-2)\ , , ,^, 
ydr"^ r dr 7'^ / 

It will be shown that quasi-exact solutions to the potential 

V{r) ^ Air + Aar^ + X^r"^ (4) 

can be obtained in similar fashion to the polynomial potentials of degree 4fc + 2 mentioned previously. It should be 
noted that the potential (4) is directly related to the quantum double well anharmonic oscillator in an external field 
studied in [29, 30]. 

In order to solve the Schrodinger equation (2) with the potential (4), we use the ansatz: 

i?,(r) =r'$(r)exp(-ar~/3rV3), (5) 

where a ^ and /? > are real parameters that arc to be determined. Similar to the polynomial potential of degree 
4fc + 2 cases, we search for possible solutions to the potential (4), of which $(r) is a finite polynomial. In this case, 
because the potential (4) is non-singular, the ansatz (5) satisfies the standard physical requirement of normalizability 
of bound states. [^^1 Namely, the ansatz (5) is finite at r = and always normalizable for r S [0, oo) even when I = 0, 
of which the situation is similar to that in quasi-exactly solvable cases of the polynomial potential of degree 4fc -t- 2, 
but different from that for singular potentials discussed in [31]. 

By substituting the ansatz (5) into (2), the differential equation for $(r) is 

-r'^ + i2Pr^ + 2ar-N-2l + l)'-^ + 
dr'^ dr 

[(2Ai + {N + 21 + l)l3)r^ - {2E + a^)r + {N + 21- l)a] $(r) = 0, (6) 

in which the real parameters a and j3 should satisfy 

/3 = V2AI, a = \2l\f2M (7) 



in order to keep the polynomial quadratic in r in front of $(r) in (6). In search for polynomial solutions of Eq. (6), 
we write 



'^{r)^^bjr\ 



(8) 



where {hj\ are the expansion coefficients to be determined. Substitution of (8) into Eq. (6) yields a matrix equation 
for the expansion coefficients {6j} and the corresponding eigenenergies E. In general, the resulting matrix equation 
is infinite dimensional. However, it can be shown that the matrix truncates to one of finite order if the following 
condition is satisfied: 



Ai = -{N + 2l + 2m + l)/3/2 for m = 0, 1, 2, 



(9) 



And in this case <I>(r) reduces to a polynomial of degree m. It follows that the matrix equation for determining the 
expansion coefficients {bj} in (8) is of the following form: 



where 



' (N + 21 - l)a 



-2ml3 




~(N + 21 - 1) 
(N + 21 + l)o 
-2E ~ a^ 
-2(m - 1)3 



Fv = 0, 



■ ■ ■ 

-2(iV + 20 ■ ■ ■ 

{N + 21 + 3)a -3{N + 2Z + 1) 

- 2E - a^ {N + 21 + 5)« 



(10) 



-4{N + 21 + 2) 



V 



-m(iV + 21 + 771-2) 

{N + 21 + 2m - l)a 

-2E - op- 



(11) 



is a matrix with tti + 2 rows and m+ 1 columns and v = (6o, hi, ■ ■ ■ , bm)"^ , in which b""" stands for transposition of b. 
With some further constraints on the parameters, solutions of (10) provide quasi-exact solutions of (6). The first row 
of F provides for two such conditions: (1) N = I — 21 or (2) N ^ 1 — 21 and 6i := ab^. The first condition is only valid 
for TV = 1 case with I = 0, while the second condition is valid for A^ > 1 cases. Therefore, when N ~ 1, the following 
eigen-equation determines possible solutions for Eq. (6): 



Pv= {2E + a^)v, 



(12) 



where 



2a 

-2m;3 

-2(r,i-l)/3 



(2™ - 2)c, 


-2/3 



(13) 



is a {m + 1) X (m + 1) quadri-diagonal matrix. Since the matrix P is non-symmetric, its eigenvalues for /3 > and 
any real a may be complex and therefore should be discarded. The remaining solutions with real eigenvalue of P are 
possible physical solutions to the problem. Therefore, the number of physical solutions in this case may be less than 
m + 1. This case is similar to the quantum double well anharmonic oscillator in an external field discussed in [29, 30]. 
However, since parity is always conserved for the potential (4), the potential (4) should be written explicitly as 



Vix) 



XlX + \2x'^ + X^x"^ , 
— XlX + X2X^ + X^x^ 



x>0 
x<0 



(14) 



when A^ = 1. It is clear that the A^ = 1 case provides part of solutions for the potential (14) on the full line with 

— OO < .T < +00. 

When TV > 1 , the eigen-equation of v becomes 



Qv = o, 



(15) 



where 



Q = 



(JV + 2! + l)a 



-2(m - 1)^ 



0- ■ ■ 

-2(JV + 21) ■ ■ ■ 

(JV + 2i+3)a -3(iV+2i + l) 

- 2S - q2 {N + 21 + 5)a 



-4(JV + 2i + 2) 



(JV + 21 + 2m - 3)q 



-4/3 



,(JV + 2i + m - 2) 
' + 21 + 2m ^ l)c, I 



(16) 



which is also a (m + 1) x [jn + 1) quadri-diagonal matrix, together with further constraint on parameters a and /3: 



/3 = -(S+-a^) 



1 2\ ^" 



Um-l 



(17) 



In this case, the condition to > 1 must be satisfied. Furthermore, because the parameter /3 depends not only on a, 
but also on eigen-energy E according to (17), similar to the N = \ case, whether A2 < or A2 > should be satisfied 
in order to keep solutions of (15) and (17) physical needs to be fixed on a case-by-case basis. In the following, we list 
some eigen-energies for small vn and the corresponding eigenvector v for iV = 1 and N >2 cases separately. 

When A^ = 1, there is only one eigenvalue with Em=o = —■^oi^, for which both A2 > and A2 < are allowed; 
there are two eigenvalues with 



p(i) 



v/-A2, 6i' 



x/^ 



"0 ' 



E. 



(2) 



4^2 +v/^, b[ 



(1) 



(,(1) 



V^A^ 



(18) 



Thus, A2 < must be satisfied for m = 1 as shown in (18). 

When N > 2, there is only one eigenvalue with E„i=i ~ ^a^(A^ + 21) with 61 ~ a6o and the constraint /? ~ 
— ^{N + 21 + l)a^, for which the condition A2 < must be satisfied; there are two independent solutions for m = 2 
with cither 



^mi2 = - o (^ + 2/ + 5 + ^/{N + 21 + 1){9N + 181 + 25))a'^, 



^(1) _ 5JV + 10^ + 5 + v/(JV + 21 + 1){9N + 181+25) 2(1) (i) _ (1) 
b2 - 8(iV + 20 " '^0 . ^1 " «^o , 



{N + 21 + 1 + ^(7V + 21 + 1){9N + 181 + 25)){5N + 10^ + 5 + y^{N + 21 + 1){9N + 181 + 25))a3 
^- 64(iV + 20 ' (^^"^ 



or 



El^U = - o (^ + 2/ + 5 - y^iN + 21 + 1){9N + 181 + 25))a^, 



(2) _ 5iV + 10Z + 5-v/(JV + 2/ + l)(9JV + 18; + 25) (2) (2) _ (2) 
2 " 8(7V + 20 « (^0 ' '^1 - "(^o ' 

(iV + 2/ + 1 - ^/{N + 21 + 1){9N + 181 + 25))i5N + 10^ + 5 - ^y{N + 21 + 1){9N + 181 + 25))c 



(1%) 



64:{N + 2l) 

Solving equations (12) or (15), we can also evaluate higher (lower) excited states with conditions that the parameters 
should satisfy. As the results arc quite complicated, they will not be provided here. 

In order to keep physical solutions only, one may project the matrices P and Q onto the physical subspace. Let 
v'*-* be i-th eigenvector with the corresponding eigenvalue to be real, and A be the projection operator, which can be 
expressed as 



A = ^vW(v«)t. (20) 



Then, instead of (12) and (15), physical solutions for A^ = 1 and A^ > 1 cases can be obtained respectively from the 
following equations: 

APAv= (2£; + a^)v, (21) 

and 

AQAv = 0. (22) 

In addition, similar to the electrostatic interpretation of the location of zeros of the Heine- Stieltjes polynomials'"^^! 
and the polynomials for the sextic anharmonic oscillator, [^"^l the polynomial $(r) satisfying Eq. (6) can also be 
determined in the same way. Let the polynomial $,„(r), up to a normalization constant, be expressed in terms of its 
zeros as 



'^.n{r) = Y[{r-r,). (23) 

Since there is no node for m — 0, and the first and second derivatives of 'J'ol^) a-re zero, 

Ai = — {N + 2l + l)/3, Eo^ — a^ and A^ = 1 - 2/ (24) 

must be satisfied according to Eq. (6), which is consistent with solution for TV = 1 and m = case reported previously. 
Since the second derivative of $i(r) is zero, the condition 

2/3r?+2a-(A^ + 2/-l)— =0 (25) 

must be satisfied. In this case, $i(r) = r — ri. It is obvious that there are two solutions of Eq. (25) with ri = ±^/^ 
when iV = 1 — 2L Once can check that there is indeed a solution of Eq. (25) with J'l = — — when N ^ 1 — 21 with 
condition f3 = —^{N + 21 + l)a^, which must be satisfied in this case for Eq. (6) with r ^ ri. For m > 2, one can 
write 

<i>"(rO V- 2 






H-^. P6) 



at any zero r^. Then, Eq. (6) at any zero r^ becomes 

i:^-2,,^-2. + ^^±ii^.O ,27, 

for i — 1, 2, • • • , r7i, of which similar equations for the generalized Lame functions arc known as the Niven 
equations. [^^1 Eq. (27) provides with important relations to determine the polynomial $m(r). Once the polyno- 
mial <^m{'r) is determined, one can get the corresponding energy and further constraints on the parameters from Eq. 
(6). The results are similar to the solutions of the sextic case^^'^l and closely related to Bethe ansatz solutions of 
Lipkin-Meshkov-Glick model!^'^'^''! and the nuclear pairing problem. ['^^'•^^1 

In summary, the 0{N) invariant quartic anharmonic oscillator is shown to be exactly solvable if the interaction pa- 
rameter satisfies some special conditions. The problem is directly related to that of a quantum double well anharmonic 
oscillator in an external field. A finite dimensional matrix equation for the problem is constructed explicitly, along 
with analytical expressions for some excited states in the system. The corresponding Niven equations for determining 
the polynomial solutions for the problem are given. It is expected that an ansatz similar to (5) may also be possible 
for finding quasi-exactly solvable cases for anharmonic oscillators with polynomials of degree 4fc for fc > 1, which add 
another quasi-exactly solvable class of anharmonic oscillators. These results should be useful to test the validity of 
various approximation method, to probe the true nature of the solution, and to reveal its asymptotic behavior. 
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